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We present a two-flavor dynamical QCD calculation of the B meson B parameters and decay constant. We use 
NRQCD for heavy quark and the nonperturbatively 0(a)-improved Wilson action for light quark at /3 = 5.2 on 
a 20^ X 48 lattice. We confirm that the sea quark effect increases the heavy-light decay constant, while estimate 
of its magnitude depends significantly on the fitting form in the chiral extrapolation. For the B parameters, on 
the other hand, we do not find a significant sea quark effect. The chiral extrapolation with logarithmic term is 
examined for both quantities and compared with the prediction of ChPT. 



1. Introduction 

In recent lattice calculations, there have been 
several indications that the presence of sea quark 
increases the B meson decay constant signifi- 
cantly, i.e. by about 10-15% [|]. Then, the study 
of the sea quark effect on the B parameters has a 
great phenomenological importance, because the 
relevant physical quantity to the B — B mixing is 
the combination /^Bb. 

In this work we perform unquenched lattice cal- 
culations of both fs and Bb using a consistent 
formulation. We use the nonperturbatively 0(a)- 
improved Wilson quark action for both sea and 
valence light quarks, and the NRQCD action for 
heavy quark. Using the same dynamical QCD 
configurations, a search for sea quark effects in 
the light hadron spectrum is also being performed 
j2j. Quenched calculations of /b and Bb using 
the NRQCD action have been presented previ- 
ously |,|. 

In full QCD, the effect of virtual pion loops 
leads to chiral logarithms, a specific dependence 
of physical quantities on the sea quark mass. 
Hence, it provides a non-trivial test of the sea 
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quark effect in lattice simulations. For /b and 
Bb chiral perturbation theory (ChPT) is exten- 
sively studied in |^,^ including its extension to 
the quenched and the partially quenched cases, 
where sea and valence quark masses can be dif- 
ferent. We test the chiral logarithm using our 
numerical data. 

2. Simulation details 

We generate a set of unquenched gauge con- 
figurations with the plaquette gauge action and 
two-flavors of nonperturbatively 0(a)-improved 
Wilson fcrmions Simulations are performed 
on a 20'^ x 48 lattice at f3— 5.2, which roughly 
corresponds to the lattice spacing a 0.1 fm. 
We take five different sea quark masses covering 
mps/my =0.6-0. 8. For comparison we also per- 
form a quenched simulation at /3=6.0 on the same 
lattice size. Further details of our simulations are 
given in 

We use the NRQCD action for heavy quark in- 
cluding all l/rriQ corrections consistently in the 
action and operators. In the calculations of /b 
and Bb we follow the methods adopted in the pre- 
vious works 1^,1). In particular, the B parameters 
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Figure 1. Chiral extrapolation of 



Figure 2. The 1/Afp dependence of 



are computed with four different prescriptions, aU 
of which have equivalent accuracy up to higher 
order contributions we neglect. We may estimate 
the systematic uncertainty taking the differences 
among them [Q. 

3. Decay constant 

Figure shows the chiral behavior of the heavy- 
light decay constant <i>j^ ~ fBy/n^B with mg ~ 
nib. Both horizontal and vertical axes are nor- 
malized with a proper power of vq and converted 
to the physical unit using the lattice result for 
rorup combined with physical input for mp. If 
we compare the unquenched result (filled trian- 
gles) with the quenched data (open diamonds), 
we clearly find the increase of about 15% with 
the unquenching as found previously |^ 

We consider the following three forms for the 
chiral extrapolation 
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J. The "log" form is motivated by 
the ChPT result for heavy- light mesons For 
two-flavor QCD, it predicts 

2 3 1 + 3ff^ 

with / the pion decay constant. The coupling g 



describes B*Btt interaction, which may be eval- 
uated from D* decays using heavy quark sym- 
metry. If we take a range g=0.2-0.7 rather con- 
servatively covering D* — > Dtt decay |^ and the 
recent D* width experiment then we obtain 
r2c2=-0.5(2) GeV-2. 

We find that both "qua" and "log" fits describe 
our unquenched data well, while the linear fit 
gives somewhat poorer than the other two. 
For the "log" fit we obtain rgc2=-0.4(3) GeV-^, 
which is consistent with the above phenomeno- 
logical estimate. 

On the other hand, in quenched QCD, no de- 
viation from linearity is observed in our data, 
and the three fit forms lead to almost identical 
result in the chiral limit. This behavior contra- 
dicts the prediction of quenched ChPT (QChPT) 
1^,^, which gives even stronger logarithmic de- 
pendence. There are, however, many unknown 
parameters in QChPT, so it is still possible to 
find a parameter region which mimics the linear 
behavior in our light quark mass regime. 

Since we are not able to decide the fit form with 
our current unquenched data, the uncertainty in 
the chiral limit is quite large. Figure |2| shows the 
1/Mp dependence of ^fs^ at the physical light 
quark mass. At the physical B meson mass, we 
observe that the two-fiavor result (filled symbols) 
is 14-25% larger than the quenched result (open 
circle), depending on the fit form. 
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Figure 3. The 1/Mp dependence of i?B^(mf,). 

4. B parameters 

We also examine the chiral behavior of Bb by 
applying the fit forms in ([^), but, because the 
light quark mass dependence itself is quite small, 
we do not observe any significant dependence on 
the fit functions. A weak dependence on sea 
quark mass is also suggested by ChPT, which 
gives 
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for Bb- 

Figure || shows the 1/Mp dependence of 
BBijnb) for both quenched (open symbols) and 
unquenched QCD (filled symbols). Four different 
symbols correspond to different prescriptions de- 
fined in [Q. Within the uncertainty indicated by 
the different prescriptions, we find no significant 
effect of unquenching. Bs^ (fTib) is also analyzed 
in a similar manner and we arrive at the same 
conclusion. 

5. Conclusions 

Our preliminary results in two-fiavor QCD are 
summarized as follows, 

fB, = 190(14)(07)(19)MeV, 

IbJIb, = 1.184(26)(20)(15), 

BBA^b) = 0.872(39)(04)(73), 

BbJBb, = 0.999(12)(04), 

BsAmt) = 0.858(33)(07)(72), 



where errors are from statistics, chiral extrapola- 
tion, and systematics. We discard the linear fit 
from the final results because of its poor value. 
The central value is obtained by taking an aver- 
age over "qua" and "log" results, and their differ- 
ence is taken as a systematic error from the chi- 
ral extrapolation. Some of the systematic errors 
in the B parameters are estimated by the scat- 
ter with different prescriptions as in Q. Other 
uncertainties are estimated by a naive dimen- 
sional counting, assuming Aqcd=500 MeV and 
a^=ay(l/a)=0.267. 

Our data for and Bb are consistent with 
the prediction of ChPT. This is different from the 
light hadron sector where we find that ChPT fails 
to reproduce the lattice data for pion mass and 
decay constant |^]. To investigate these issues 
we are currently accumulating further statistics, 
which also helps to reduce the statistical error in 
our results. 
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